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ABSTRACT.  7  A.  systematic  development  from  the  Dyson  equation  and  the 
two-frequency  Bethe-Saltpeter  equation  of  a  two-frequency  radiative  trans¬ 
fer  equation  suitable  for  pulsed  waves  in  the  presence  of  a  random  distri¬ 
bution  of  absorptive  discrete  scatterers  with  pair  correlations  is 


presented.  The  main  strength  of  the  radiative  transfer  theory  expounded 
here  stems  from  the  fact  that  it  is  applicable  under  conditions  of  large- 
angle  scattering,  statistical  inhomogeneities  and  statistical  anisotropies. 
It  accounts,  also,  for  regular  refraction  (variable  scatterer  density,) 


absorption  and  frequency  offsets. 


I.  INTRODUCTION.  Multiple  scattering  by  a  random  distribution  of 
discrete  scatterers  has  been  studied  extensively  over  the  past  thirty  five 
years,  primarily  because  of  its  relevance  to  a  large  number  of  pressing 
applied  problems  that  arise  in  radio  physics  and  engineering.  Fundamental 
work  on  multiple  scattering  of  scalar  waves  by  a  distribution  of  uncorrelated 
scatterers  was  initiated  by  Foldy  [1]  and  has  been  developed  further  by 
Lax  [2]  and  Twersky  [3].  Based  on  these  original  contributions,  a  great 
number  of  applications  have  appeared  in  the  literature  [cf.,  for  example. 
Refs.  4-7]  involving  both  the  coherent  field  and  the  incoherent  intensity. 

The  problem  of  assessing  the  multiple  scattering  effects  on  scalar  and 
vector  waves  in  the  presence  of  random  distributions  of  correlated 
scatterers  is  more  challenging.  Basic  contributions  along  this  direction 
have  been  made  by  Twersky  [8],  Bringi  et  al.  [9]  and  Tsang  and  Kong  [10]  in 
connection  with  the  coherent  field,  and  by  Barabanenkov  [11],  Barabanenkov 
and  Finkel'berg  [12]  and  Watson  et  al.  [13]  in  connection  with  the  incoherent 
field  Intensity. 


The  motivation  for  our  work  is  based  on  the  absence  of  a  second-order 
statistical  theory  for  studying  pulse  propagation  through  random  distributions 
of  correlated  scatterers.  Such  a  theory  is  necessary  for  the  development 
of  predictive  models  pertinent  to  (1)  pulsed  electromagnetic  propagation 
through  complex  natural  media  (e.g.,  rain,  fog,  sandstorms,  vegetation); 

(2)  the  analysis  of  obscuration  and  detection  techniques;  (3)  remote 
sensing  and  identification  of  aerosol  clouds;  (4)  scattering  clouds 
consisting  of  complicated  individual  scatterers  (lossy,  anisotropic, 
frequency-sensitive,  possibly  aligned,  pair-correlated)  with  gross  macro¬ 
scopic  structure  (finite  extent  and  variable  number  density). 
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Our  specific  goal  in  this  exposition  is  to  derive  systematically 
a  radiative  transfer  equation  for  the  two-frequency  incoherent  intensity 
function  (a  quantity  related  to  the  two-frequency  mutual  coherence  function 
and*  hence*  to  second-order  pulse  statistics.)  The  derivation  is  limited 
to  scalar  pulsed  waves  propagating  through  a  random  distribution  of 
absorptive  scatterers  with  pair  correlation}. 

The  proposed  radiation  transport  equation  is  based  on  the  Dyson  and 
Bethe-Saltpeter  equations  at  the  level  of  the  direct  interaction  and 
ladder  approximations,  respectively.  If,  ii  addition  to  pair  correlations, 
the  assumptions  are  made  that  the  number  of  scatterers  is  large  and  the 
average  distance  between  any  two  scatterers  is  large  compared  with  a 
reference  wavelength,  the  Dyson  and  Bethe-Saltpeter  equations  are  analogous 
to  those  associated  with  a  continuous  random  medium  with  fluctuations  of 
the  permittivity  which  are  distributed  according  to  a  normal  law  and  with 
a  deterministic  profile  directly  linked  to  the  number  of  scatterers  per 
unit  volume. 

The  transition  from  the  Dyson  and  Bethe-Saltpeter  equations  to  the 
two-frequency  radiative  transfer  equation  is  effected  by  a  continuous 
stochastic  transport  theory  that  was  originally  introduced  by  Barabanenkov 
at  al.  [14]  and  subsequently  extended  to  the  two-frequency  context  by 
Besleris  and  Kohler  [15-17].  As  in  their  case ,  our  derivation  differs 
markedly  from  the  usual  procedures  for  obtaining  classical  radiation  transport 
equations;  the  latter  rely  mostly  on  consic erations  of  energy  balance,  with 
no  explicit  "microscopic"  interpretation  given  to  the  extinction  and 
scattering  coefficients. 

In  the  following  we  present  a  sketch  of  the  proposed  two-frequency 
radiation  transport  theory,  with  primary  eiiphasls  on  the  underlying  assump¬ 
tions  and  the  physical  interpretation  of  tl  e  various  terms  entering  into  the 
final  transport  equation. 

II.  FIRST  AND  SECOND  ORDER  COHERENCE  FUNCTIONS  FOR  A  RANDOM  DISTRIBUTION 
Of  PAIR-CORRELATED  SCATTERERS.  PARI  A:  ANISOTROPIC  SCATTERING.  The  deriva¬ 
tion  of  the  Dyson  equation  for  the  mean  field  and  the  Bethe-Saltpeter 
equation  for  the  mutual  coherence  tensor  under  conditions  of  anisotropic 
scattering  and  in  the  presence  of  pair  correlations  among  scatterers  is 
based  on  the  Twersky  procedure.  The  basic  underlying  assumptions  are  the 
following:  (1)  We  ignore  third-order  scattering  by  two  scatterers,  fourth- 
order  scattering  by  three  scatterers,  etc.  (essentially  the  Twersky 
assumption);  (2)  All  scatterers  have  the  same  shape,  size  and  orientation 
distributions;  (3)  We  consider  only  pair  correlations  and  neglect  all 
contributions  from  higher-order  correlations;  (4)  The  number  of  scatterers 
in  a  volume  is  infinite. 


We  resort,  also,  to  the  following  notational  definitions:  (1)  £a: 
incident  electric  field  at  position  r  ;  (2)  E3 :  total  electric  at  r  ; 
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(3)  g  a  :  electric  field  at  £  caused  by  the  jth  particle;  (4)  <<•>>: 

j  a 

configurational  averaging  (over  size,  shape  and  orientation  distributions) ; 
(5)  p(r)j  particle  density  function;  (6)  B'(r  ,  r.) :  pair  correlation 
function. 

Under  the  aforementioned  assumptions,  the  Dyson  equation  for  the 
coherent  vector-valued  electric  field  assumes  the  form 


<E*>  -  +  /  dr^  «gjS»  <E^>  P(£j) 

+  /  d£j  /  **  <<8j3>>  «Gkj»  <Ek>  B’O^,  rfc)  ; 


(2.1a) 


«Gka»  =  «gka»  +  /  dr^  «g 


,a>>  <<Gfc£>>  p<££> 


+  /  dr  /  dr  «g  a»  «G  A»  «G  “»  B’(_r.,  r  )  .  (2.1b) 


On  the  other  hand,  the  Bethe-Saltpeter  equation  for  the  mutual  coherence 
tensor  becomes 

<EaEb*>  *  <Ea>  <Eb*>  +  /  d_Tj  «G^a»  <E^  E^*>  «G^b*»  p(r^) 

+  /  dr,  /  dr  «G.a»  <E*  E®*>  «G  b*»  B'(r,,  r  )  .  (2.2) 

A  detailed  derivation  of  Eqs.  (2.1)  and  (2.2)  is  given  in  Ref.  18. 

The  former  is  at  the  level  of  the  direct  interaction  approximation,  whereas 
the  latter  is  at  the  level  of  the  ladder  approximation.  Both  are  analogous 
to  equations  associated  with  a  continuous  random  medium  with  fluctuations 
of  the  permittivity  which  are  distributed  according  to  a  normal  law  and  with 
a  deterministic  refractive  profile  directly  linked  to  the  number  of 
scatterers  per  unit  volume.  It  should  be  noted,  however,  that  no  counter¬ 
part  to  the  second  (collapsed)  term  on  the  right-hand  side  of  (2.2)  exists 
in  the  continuous  random  medium  case.  In  the  absence  of  pair  correlations, 
i.e. ,  B'(r, ,  r  )  *  0,  (2.1)  and  (2.2)  reduce  to  the  equations  derived 

“j  —in 

previously  by  Twersky  [cf.  Ref.  19]. 

The  vector-valued  Dyson  equation  (2.1)  and  the  tensor-valued  Bethe- 
Saltpeter  equation  (2.2)  are  the  basic  equations  for  deriving  a  tensor¬ 
valued  radiative  transfer  equation  for  vector  waves.  Such  a  derivation 
is  very  complicated  and  will  not  be  undertaken  in  this  paper.  Instead, 

Eqs.  (2.1)  and  (2.2)  will  be  "scalarlzed"  in  the  next  section.  (A  set  of 
assumptions  sufficient  for  such  an  approximation  i*  given  in  the  Appendix.) 
The  resulting  equations  will  form  the  basis  for  deriving  a  scalar  radiation 
transport  theory  in  Sec.  IV. 


Trim 


III.  FIRST  AMD  SECOND  ORDER  COHERENCE  FUNCTIONS  FOR  A  RANDOM 


DISTRIBUTION  OF  PAIR- CORRELATED  SCATTERERS ■  PART  B:  ISOTROPIC  SCATTERING. 


In  the  following,  we  shall  assume  that  the  scattering  channel  Is  tenuous, 
l.e.,  the  distance  between  any  two  scatterers  Is  much  greater  than  a 
reference  wavelength.  He  shall  also  assume  scalar  Isotropic  scattering. 

In  the  electromagnetic  case  this  approximation  arises  if  individual 
scatterer  dimensions  are  small  compared  to  wavelength  and  if  "gross" 
depolarisation  effects  are  neglected.  In  the  case  of  acoustic  wave 
propagation,  isotropic  scattering  takes  place  when  individual  scatterer 
dimensions  are  small  compared  to  wavelength. 

Under  these  assumptions,  one  deals,  essentially,  with  a  scalar  wave 
theory.  The  Dyson  equation  (2.1)  simplifies  considerably  and  is  rewritten 
below  in  a  form  suitable  for  our  work  in  Sec.  IV: 


(V2  +  k2  +  4wf  p(r)]  <E(r,  k)> 

-  -4irf  /  dr'  «G'  (r ,  r ' ,  k)»  <E(r',  k)>  B'(r,  _r’)  ; 


(3.1a) 


[V2  +  k2  +  4irf  p(r)]  «G'(r,  r ' ,  k)» 

-  -4irf  «(jr  -  ir')  -  4irf  /  dr"  «G'  (jr,  r",  k)»  «G’  (r",  r' ,  k)»  B'(r, 

(3.1b) 


The  notation  is  identical  to  the  one  used  in  the  previous  section,  except 
that  E  and  Gf  are  now  scalar-valued.  The  quantities  k  and  f  are  respectively 
the  wavenumber  and  a  configurationally  averaged  scalar-valued  scattering 
coefficient.  The  latter  is  independent  of  _r;  however,  it  may  depend  on  k 
and,  in  general,  is  complex  by  virtue  of  the  absorptive  properties  of  the 
scatterers. 

The  scalar-valued  Bethe-Saltpeter  equation  for  the  two-frequency 
mutual  coherence  function  r' (r^,  r^*  ^  H  <E(£1»  kj)E*(r2,  k2^  can  be 

written  as  follows:  1 


1<7r  2“Vr22)  +<kl2  “  k22)  +  4irf  PiLl)  ‘  4irf*  p(i2n  T'^V  ^2’  kl»  k2} 
•  -4w{p(r^)f  «G'*(r2,  r1#  k2)»  r'^,  rx,  k2) 

“P(E2>f*  <<G,^l»  -2*  klJ>>  r'<l2*  -2*  kl»  k2)J 


-4*  /  dr’  [f  B'Q^,  r’)  -  f*  B’  (r* ,  r2)H«C,CE1.r,>  k1>»  r’(r’,  r. v  k1#  k 


+  «C'*(r2,  £',  k2)»  r’^,  r',  k^  kj)] 


(3.2) 


The  derivation  of  (3.2)  [cf.  Ref.  18]  Incorporates  a  narrowband  pulse 
assumption.  In  this  case,  the  scattering  coefficient  f  depends  only  on  the 
prescribed  carrier  frequency. 

The  functional  forms  of  Eqs.  (3.1)  and  (3.2)  are  analogous  to  those 
derived  by  Beaieris  and  Kohler  [15,16]  for  acoustic  wave  propagation  in  a 
continuous  random  medium  characterized  by  Gaussian  fluctuations.  As  in  the 
case  of  the  more  general  equation  (2.2),  no  counterpart  of  the  first 
(collapsed)  term  on  the  right-hand  side  of  (3.2)  exists  in  the  continuous 
random  case.  For  k^  *  k2  and  in  the  absence  of  pair  correlations,  i.e., 
B'(r2,  jr2)  -  0,  Eqs.  (3.1)  and  (3.2)  reduce  to  relationships  already 

available  in  the  literature  [cf.  Ref.  19]. 

.  TWO- FREQUENCY  RADIATIVE  TRANSFER  EQUATION  FOR  A  RANDOM  DISTRIBUTION 


OF  ABSORPTIVE  SCATTERERS  WITH  PAIR  CORRELATIONS:  SCALAR  WAVE  THEORY.  In  the 


Bethe-Saltpeter  equation  (3.2)  for  r  (r ^ ,  £2»  ki»  w®  intr°duce  center-of 

mass  and  difference  coordinates  and  wavenumbers,  viz.,  R  =  (r^  +  £2)/2, 

r  5  _r2  -  £2;  kg  =  (k^  +  k2)/2,  k^  =  k^  -  k2,  and  use  the  notational 

definitions  f'fej,  £2,  kj,  k2)  =  T(R,  r,  kg,  kd),  B' (r^,  £2)  =  B(R,  r), 

«C,(r1,  £2,  k)»  =  «G^R,  £,  k)»,  and  M(R,  £,  k)  ~  f«G(R,  £,  k)»  B(R,  £) 

We  introduce,  also,  the  Fourier  transform  pairs  r(R,  £,  kg,  k^)  ■* — *•  f(R,  k, 

kg,  kd>»  B<*»  I)  •(£.  «>»  aad  **(£»  £.  k>  M(R,  £,  k)  -  M'(R,  k,  k)  + 

i  M"(R,  k,  k).  It  should  be  noted  that  M(R,  £,  k)  is  analogous  to  the  "mass 
operator"  entering  into  the  Dyson  equation  in  the  case  of  smoothly  inho¬ 
mogeneous  media.  Furthermore,  the  quantity  f (R,  jc,  kg,  k^)  is  the  two- 

frequency  extension  to  the  phase-space  Wigner  distribution  function. 

We  consider  next  smoothly  inhomogeneous  media  for  which  T(R,  jr,  kg,  k^) 

B(R,  r)  aad  «G(R,  £,  k)»  vary  slowly  with  respect  to  the  sum  variable  R, 
and  rapidly  with  respect  to  the  difference  variable  £.  We  make  also  two- 
further  assumptions:  (1)  the  ratio  of  difference  to  sum  wavenumbers  is 
small  cosq>ared  to  unity,  i.e.,  jk^/k  |  «  1;  (2)  the  scattering  and  regular 


losses  are  small  but  not  negligible.  Within  the  fraawwork  of  these  restric¬ 
tions  scattering  becomes  significant  [cf.  Ref.  18  for  details]  on  the 
"energy"  surface 


*'  (*♦  J£»  kg)  2  y[ic2  -  kg2  -  4irfR  p(R)  -  4ir  M'  (R,  jc,  kg)  ]  - 


(4.1) 


which  is  Independent  of  regular  and  scattering  losses.  The  quantity  in 
(4.1)  denotes  the  real  part  of  the  complex  scattering  coefficient  f. 

In  the  general  case  of  statistically  anisotropic  pair  correlations  of 
the  scatterers  we  seek  a  solution  in  the  form 


f(R,  jc,  k,  k.)  -  f  (R,  k,  k  ,  k.)  +  f .  (R,  k,  k,  k.)  .  (4.2) 

The  first  term  on  the  right-hand  side  of  (4.2)  is  the  coherent  part  of  the 
two-frequency  Wigner  distribution  and  is  directly  related  to  the  solution 
of  the  Dyson  equation  (3.1).  The  incoherent  part  of  the  Wigner  density 
function,  on  the  other  hand,  is  chosen  as  follows: 


£»  k8.  V  -  kg|7K  H’(R,  £,  kg)||VKH'(R,  £,  kg)|“3 

*  fi[H*  (R,  £,  kg)]  I(R,  £,  kg,  kd)  ;  s  -  £/ic  .  (4.3) 


The  quantity  I(R,  s,  kg,  kd)  is  the  two-frequency  Incoherent  "ray"  intensity 
at  the  point  R  and  in  the  direction  of  the  unit  vector  js. 


Let,  next,  kgf^  (R,  t,  kg)  denote  the  value  of  k  for  which  H'(R,  k£. 


kg)  ■  0,  and  define  an  effective  index  of  refraction  as  follows: 


neff (5»  £»  k8>  *  lvK  H'CR.  «£»  k8)!/k8  ;  x  “  keff(*»  £»  k8>  •  (4*4) 


Let,  finally,  0(R,  s,  k#)  be  the  angle  between  the  direction  of  the  group 

and  phase  velocities.  With  these  definitions  in  mind,  the  incoherent  ray 
intensity  I(R,  ja,  kg,  k^)  is  found  to  obey  the  two-frequency  radiative 

transfer  equation 


®8ff(R»  £»  k8)  -£[  {I(R»  £.  kg»  kd)l«>»  9(R»  £»  k8)|_1  neff ’  -*  ks» 


-  (-4  fj  p(R)  -  4w  M"[R,  kgff  (R,  s,  kg)s,  R.l+iJ  kgkj} 


x  k81  n8}f  <*»  £»  k8>  I<1*  £»  kd> 


+  2 1 f | 2  kg2  /  ds*  kjff  (R,  s’,  kg)ngff(R,  £,  k8>neff&»  *.'»  k8> 


*  | cos  0(R,  s' ,  kg) |-1{$[R,  kgff(R,  £,  kg)js  -  kgff(R,  £,  kg)^' ]  +  p(R)} 

x  I(R,  s’,  k  .  k.) 

—  —  s  d 

+  2 1  f  | 2  ngff(R,  £,  kg)  |  cos  0(R,  £,  kg)  |  /  djc’  {<t>[R  -  kgff(R,  s,  kg)£  -  jc1 '  ] 

R3 

+p(R)}  fQ(R,  £»,  kg,  kd)  ,  (4.5) 


where  ft  denotes  the  range  of  over  the  surface  of  a  unit  sphere. 

Equation  (4.5)  for  the  two-frequency  incoherent  intensity  I(R,  s,  k  ,  k.) 

—  —  s  d 

is  the  main  result  of  this  paper.  Zn  interpreting  this  equation,  we  should 
note  the  following:  The  left-hand  side  of  the  equation  is  a  convective 
term;  the  ray  paths  correspond  to  an  effective  medium  determined  by  the 
density  of  the  scatterers  p(R),  the  scattering  coefficient  f  and  the 
spatial  correlation  function  of  the  scatterers.  [The  quantity  dl  denotes 
the  differential  of  a  curvilinear  ray  passing  through  the  point  R  in  the 
direction  V  H' (R,  jc,  k  ).]  The  first  and  second  terms  on  the  right-hand 

side  of  (4.5)  are  due  respectively  to  regular  and  scattering  losses;  the 
third  term  arises  because  of  frequency  offsets;  the  fourth  one  is  the 
scattering  term;  finally,  the  last  one  is  the  source  term,  representing 
the  "feeding"  of  the  incoherent  intensity  by  the  coherent  part  of  the  Wigner 
distribution  function. 

In  order  to  compute  the  second-order  pulse  moment  <E(r^»  t^E^r^*  t2^> 

at  the  receiver  site,  we  must  first  find  the  two-frequency  mutual  coherence 
function  E(_rn ,  gj^)E*(£2»  f^)  and  then  perform  a  two-dimensional  Fourier 
transform  with  respect  to  and  u>2.  In  the  ladder  approximation,  the 
two-frequency  mutual  coherence  function  obeys  the  Bethe-Saltpeter  equation 

(3.2) .  The  phase-space  analog  of  the  two-frequency  mutual  coherence  function 
is  the  Wigner  distribution  f(R,  £,  k  ,  k,).  The  latter  is  decomposed  in 

(4.2)  into  a  coherent  part,  which  is  directly  linked  to  the  solution  of  the 

Dyson  equation  (3.1),  and  an  incoherent  part  directly  expressible  [cf. 

Eq.  (4.3)]  in  terms  of  the  two-frequency  ray  intensity  I(R,  s,  k  ,  k,). 

—  —  s  a 

The  latter  obeys  the  radiative  transfer  equation  (4.5). 

Not  all  the  successive  steps  outlined  in  the  previous  paragraph  need 
be  followed  for  obtaining  information  about  second-order  pulse  statistics. 

If,  for  example,  such  information  is  restricted  to  the  Incoherent  part 
of  <E(r^,  t^)E*(£2»  t 2 )>»  direct  usage  can  be  made  of  the  two-frequency 

ray  intensity  I(R,  £,  kg,  kd),  as  explained  below. 
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Let  ri(R,  £,  kg,  k^)  — ►  f^(R,  £,  kg,  k^)  denote  the  incoherent  part 
of  r(R,  jc,  kg,  k^).  Then, 

rt(l»  £,  kg,  kd)  -  /  da  kg-2  n”|f  (R,  s,  kg)k2ff(R,  s,  kg) 

X  | cos  0(R,  e,  kg) I”1  I(R,  a,  kg,  kd)exp{i  heff(R.  £,  kg)£  •  r>  ,  (4.6) 

which  establishes  a  useful  connection  between  the  photometric  ray  intensity 
I(R,  s,  kg,  kd)  and  the  incoherent  part  of  the  two-frequency  mutual  coherence 

function. 

V.  CONCLUDING  REMARKS.  A  radiative  transfer  equation  for  pulsed  scalar 
waves  in  a  random  distribution  of  pair  correlated  absorptive  scatterers 
has  been  derived  systematically  from  the  Dyson  and  Bethe-Saltpeter  equations. 
Detailed  solutions  —  both  analytical  and  numerical  —  are  presently  under 
consideration.  An  important  question  in  this  context  is  whether  controlled 
experiments  could  be  carried  out  so  that  comparisons  would  be  made  with 
theoretical  results. 

An  open  research  area  in  propagation  through  random  distributions  of 
scatterers  is  the  systematic  derivation  of  a  tensor  transport  theory  from 
the  Dyson  equation  (2.1)  and  the  Bethe-Saltpeter  equation  (2.2).  Such  a 
theory  is  of  paramount  importance  for  physical  situations  where  anisotropic 
scattering  and  depolarization  effects  cannot  be  neglected. 

AOMOWLIDGMBITS.  Research  supported  in  part  by  the  Army  Research  Office 
under  contract  Ho.  DAAG  29-79-C-0085  and  the  Office  of  Haval  Research  under 
contract  Ho.  H00014-76-C-0056 . 

tframil.  Consider  the  case  of  Born  approximation  valid  for  single 
particle  scattering  together  with  the  simplification 


a  a 


/  U-fe')  ”  U  e*p{ik(i  -  o )  •  r’}  =  /  dr' le  (r')  -  1] 
V.  Vs 


(A-l) 


In  this  expression,  V  denotes  the  volume  of  the  scatterer,  e  (r)  is  the 

8  a  a  r  ,ir 

relative  permittivity  of  the  scatterer,  and  1,  o  are  respectively  unit 
vectors  along  the  incident  and  scattered  directions.  In  this  case,  the 

tensor-valued  quantity  <<g^a>>  obeys  the  equation 


[Vt  x  Vt  x  (.)  -  k2I]  «gJ*>>  -  k2  5  6(1^  -  r^)I  , 


(A-2) 
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where 
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